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Figura 1: Parque del Milenio, Chicago. F.O. Gehry, 2004 


Las curvas del diseño asistido por ordenador (CAD, Computer Aided Design), son 
las curvas NURBS (de Non uniform rational B-spline). Son las curvas que se utilizan 
para crear formas libres y evolucionaron a partir de las curvas de Bézier (De Casteljau 
1959, Bézier 1962) y de las curvas B-spline (Basis Spline). 

Las curvas NURBS están definidas a partir de un conjunto de puntos de control. 
Hay dos técnicas diferentes para crear curvas a partir de puntos de control. 


= Interpolación: Dado un conjunto de puntos ordenado, se genera una curva suave 
que pasa por todos esos puntos. 


= Aproximación: Se traza un polígono de control que une los puntos ordenados y 
se busca una curva que aproxime la forma indicada por el polígono de control. 


Las curvas NURBS se obtienen mediante un proceso de aproximación. 

Estos apuntes están especialmente dirigidos a alumnos del Taller Experimental Il: 
Proyectos de Arquitectura Paramétrica, y alumnos del Grado en Fundamentos de la 
Arquitectura, interesados en profundizar en las herramientas matemáticas utilizadas 
para crear curvas de forma libre (en Grasshopper y Rhinosceros). En la medida de lo 
posible, se ha intentado suavizar el rigor matemático, utilizando sólo el nivel necesario 
para trasmitir las posibilidades que ofrecen las curvas NURBS al diseñador. Se ha 
puesto énfasis en conectar el uso de Grasshopper con las matemáticas que soportan 
este programa de ordenador. 

Grasshopper es un componente de Rhinosceros 3D que se ha popularizado por su 
capacidad para generar diseño paramétrico. Es un editor gráfico de algoritmos que no 
requiere conocimientos de programación, pero permite al diseñador la construcción 
de generadores de formas tanto simples como asombrosas. 

Este es el primer volumen de la serie: 


= FORMAS LIBRES lI: Curvas NURBS. 
= FORMAS LIBRES Il: Superficies NURBS. 


1 Curvas de Bézier 


Para definir una curva de Bézier partimos de n +1 puntos Bo,... B,, planos (en 
R2) o espaciales (en I¡R*) que llamaremos vértices o puntos de control. El polígono de 
control es el polígono que obtiene tras unir el punto B; con el B;y1,1=0,...,n—1 
mediante un segmento, que podemos parametrizar por 


Bt) = (1-0 Bi +tB;, t € [0, 1]. 
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Ejemplos 1.0.1. 1. El polígono de control determinado por los puntos By = (1,1), 
B, = (3,5), B2 = (5,3), B3 = (7,4) es el de la figura 1.1. Las parametrizaciones 
de los segmentos son 


B¿(t) = (142t, 14+4t), Bl(t) = (3+2t,5-2t), Bl(t) = (5+2t,3+u), t € [0, 1]. 


2. Para puntos By = (1,1,0), B, = (1,3,5), B2= (5,-1,3), B3 = (3,4,1) en R* 
las parametrizaciones de los segmentos del polígono de control (en la figura 1.1) 
son: 

BJ(t) = (1, -2t+3,5-—5€), Bl(t) = (-4t+5,4t — 1,2 +3), 

Bi(t) = (2t +3, —5t +4,2t+1), t e [0, 1]. 


1.1. Algoritmo de De Casteljau 


Una forma sencilla de comprender cómo se obtiene la parametrización de una 
curva de Bézier, a partir de los vértices del polígono de control, es el algoritmo de De 
Casteljau, para 4 vértices o puntos de control Bo, B,, Ba, B3. Sea t un parámetro en 
el intervalo (0, 1]. Los pasos del algoritmo son los siguientes: 


1. Interpolar linealmente pares de puntos consecutivos 
Bi(t) = (1—1)Bo+tB,, 
B¡(t) = (1-t)B, +tB», 
B¿(t) = (1—1)B2+tBx, 


obteniendo la parametrización de las rectas del polígono de control. 
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Figura 1.1: Polígono de control determinado por 4 puntos en R? (izq) y R* (dcha) 


2. Interpolar linealmente pares de puntos genéricos consecutivos Bj(t), Bl(t), 
BL(t), para un mismo valor del parámetro t 


By(t) = (1 —t)B5(1) +tB,(0, 
B¡(t) = (1 —1)B¡(t) +tBa(t), 


obteniendo parametrizaciones de curvas de grado dos. 


3. Interpolar linealmente los puntos genéricos de las curvas de grado 2, para un 
mismo valor del parámetro t 


P(t) = (1 —t)By(t) +tB¡(0), 
= (1 —tBo +3(1 —t) Br +3(1 — t)t*Ba +1*Ba. 


obteniendo la parametrización P(t) de la curva de Bézier, cuyas componentes 
son polinomios de grado 3. 


Ejemplo 1.1.1. Calculamos con Maple, siguiendo este algoritmo, la parametrización 
de una curva de Bezier de grado 3 a partir de cuatro puntos By = (0,0), B, = 
(1,3), B2 = (5,4), B3 = (7, —1). Calculamos la parametrización de las rectas del 
polígono de control, B¿(t), BL(t), BL(t) que podemos ver como tres puntos genéricos 
(parametrizaciones de curvas de grado 1). Unimos estos tres puntos mediantes rectas 
B¿(t), Bi(t) obteniendo dos nuevos puntos genéricos (parametrizaciones de curvas 
de grado 2) y finalmente, la curva de Bézier viene dada por la recta que une dichos 
puntos, una parametrización dada por polinomios de grado 3 


P(t) = (3t + 98 — 5t*%, 9% — 61? — 4t*) . 


Representamos la curva junto a su polígono de control, en la figura 1.2. 
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Figura 1.2: Polígono de control de una curva de Bézier de grado 3 
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Figura 1.3: Curvas de Bézier de grado 3 


Conocidas las propiedades básicas de las curvas de Bézier un usuario puede pre- 
decir (con un poco de práctica) la forma de una curva generada por un polígono de 


control dado (véase la figura 1.3). 


1.2. Definición y propiedades 


La parametrización que proporciona el algoritmo de De Casteljau se puede ex- 


presar mediante una fórmula en la base de los polinomios de Bernstein. Dados n + 1 
., B,, definimos una curva de Bézier mediante una para- 


puntos de control By, B;,.. 


metrización, dada por polinomios de grado n 


P0)=>» Bda, DEbEl 
1=0 


con polinomio de Bernstein .,, ¿(t) = ( E ) A 


La primera propiedad, fácilmente deducible de la parametrización, es que la curva 
de Bézier contiene a los puntos inicial y final del polígono de control 


Bo=P(0) y B,=P(1). 


1.2.1. Vectores tangentes y derivadas 


El vector tangente a la curva de Bézier en el punto P(t) viene dado por su primera 
derivada P"(t) [GAS] 


Pl) = Y B,J,.L0) 


Comprobamos fácilmente que los vectores tangentes en los puntos inicial y final de 
la curva de Bézier tienen la misma dirección que el primer y último segmento del 
polígono de control, respectivamente (véase la figura 1.3). 
P*(0) =n(B1 — Bo), 
P(D)=n(B,=Ba1). 
Por otra parte, la segunda deriva P”(t) está relacionada con la curvatura [GAS]. En 
los puntos inicial y final de la curva, la curvatura depende de los tres puntos de control 
más cercanos (de las dos aristas más cercanas del polígono de control). 
P"(0) = n(n + D(Bo => 2B: + Ba), 
PU = ní(n = D(B, => 2Bn-1 + Bn-2). 
Ejemplo 1.2.1. Volvamos al caso de un curva de Bézier con 4 puntos de control. 
Derivando la parametrización obtenida en la Sección 1.1 obtenemos una interesante 
relación 
Pt) ==3(1=0*B0 + 3210) +(1-039B1 + 3(=4S+2(1—=¿))B2 + 314Ba 
=>» =3(BÍ(t) — By(t)). 
Es decir, el vector tangente a la curva de Bézier en el punto P(t) es proporcional al 
segmento Bi(t) — Bi(t) 
Si la curva de Bézier de grado 3 es espacial, los tres primeros puntos de control 


Bo, B,,B3 determinan el plano osculador en By y los tres últimos B;, Ba, B3 deter- 
minan el plano osculador en B3. 


1.2.2. Envolvente convexa 


Intuitivamente hablando la envolvente convexa de un conjunto de puntos en el 
plano se obtiene como sigue. Colocamos alfileres sobre dichos puntos y soltamos una 
goma a su alrededor, el polígono que forma la goma es la envolvente convexa. Con 
puntos espaciales podemos pensar en envolverlos con una membrana elástica. En el 
caso plano obtenemos un polígono convexo y en el espacial un poliedro convexo. En 
ambos casos, ser convexo significa que el segmento que une dos puntos cualesquiera 
está contenido en la envolvente. 

Una curva de Bézier está contenida en la envolvente convexa de sus puntos de 
control. Esto se deduce de su construcción mediante el Algoritmo de De Casteljau, que 
realiza iteradamente interpolación lineal. Este hecho aporta control sobre el diseño. 
Observamos además, que las curvas de Bézier tienen precisión lineal, en el sentido de 
que si los puntos de control están alineados, la curva de Bézier está contenida en la 
recta que determinan. Véase el ejemplo de la figura 1.5. 


1.2.3. Otras propiedades y limitaciones 


Las curvas de Bézier son invariantes por transformaciones afines. Esto significa 
que la imagen mediante una transformación afín de una curva de Bézier coincide con 
la curva de Bézier de los transformados de los puntos de control. Ejemplos de trans- 
formaciones afines son las isometrías (traslaciones, rotaciones, simetrías), semejanzas 
(composición de isometrías con homotecias) y proyecciones (véanse [CLI] y [PAHK]). 

Aunque sirven para modelar un amplio conjunto de curvas, las curvas de Bézier 
presentan dos limitaciones importantes: 


1. Un número de puntos de control elevado produce curvas que se asemejan menos 
al polígono de control, el diseñador pierde la intuición de la curva que va a 
producir. 


2. Cada punto de control tiene un efecto global sobre la curva de Bézier obtenida, 
ya que en la definición está multiplicado por un polinomio de Bernstein, con 
dominio en todo el intervalo [0,1]. Por lo tanto, si desplazamos un punto o 
incluimos un nuevo punto, cambia la forma de la curva entera. 


1.2.4. Componentes Grasshopper 


Esta primera aproximación a las curvas NURBS permite al alumno empezar a 
comprender el funcionamiento real de algunos de los componentes más importantes 
de Grasshopper. Dada una lista P de puntos de control, el componente PLine dibuja 
el polígono de control que determinan. Para dibujar la curva de Bézier utilizamos 


Figura 1.4: Componentes PLine, Nurbs, Hull y CP de Grasshopper 


Figura 1.5: (Dcha) Polígono de control de una curva de Bézier. (Izq) Envolvente 
convexa de los puntos de control de una curva de Bézier plana 


el componente Nurbs con grado D igual al número de puntos de control menos uno 
(véase la figura 1.4). En general Nurbs dibuja curvas B-spline como veremos más 
adelante. 

El componente Hull de Grasshopper devuelve la envolvente convexa de un con- 
junto de puntos del plano. En la figura 1.5 se puede ver una curva de Bézier de grado 
4, creada con Nurbs. 

Dada una curva cualquiera definida en Grasshopper, podemos visualizar sus pun- 
tos de control utilizando CP. Este componente devuelve además la lista de pesos W 
y el vector de nudos K, como veremos posteriormente. 
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2 Curvas B-spline 


Una curva definida a partir de un polígono de control depende de la estrategia de 
interpolación o aproximación que se utilice para construir la curva, depende de la base 
de funciones que se utilice para este objetivo. Utilizando los polinomios de Bernstein 
aparecen dos limitaciones importantes, el grado de la curva obtenida es siempre el 
número de vértices de control menos uno y la modificación de un vértice afecta a 
toda la curva, ya que las funciones J,,;(t) afectan globalmente a toda la curva. La 
base B-spline (de Basis spline) corrige estas limitaciones y es la que define las curvas 
B-spline. 

Aunque el arquitecto o diseñador normalmente no necesita manejar la fórmula 
que define las curvas B-spline, tener una idea del papel que juegan los ingredientes 
necesarios para definirla puede ayudar a entender su funcionamiento. Más interesante 
que la fórmula, son las distintas formas de las curvas B-spline, que resultan de las 
funciones de la base B-spline, y que varían atendiendo a la elección de un vector de 
nudos y del grado k (en la literatura también se habla de orden igual a k + 1), véase 
el ejemplo de la figura 2.1. Para definir una curva B-spline necesitamos: 


1. Los vértices del polígono de control Bg,..., B,. 
2. El grado k (orden k + 1) menor o igual que n. 


3. Un vector creciente de nudos X = (a =21,...,Un+k+2 =D), del que hablaremos 
en la Sección 2.1. 
La parametrización de la curva B-spline está definida por [Rogers] 
n+1 


P(t) = > BiaNil0), (07 < t < b 
i=1 


con base de funciones B-spline Ni ¡,...,Nn+11, 1 = 1,...,k +1 definidas de forma 
recursiva por las fórmulas de Cox-De Boor, 


l sizí¿<t< Zi 
ÑN, e A == 29 do 
at) Í O en otro caso, 
t— IN it á —- UN; t 
Nara) =é ri) Nenlt) i ESAS a, +0) SO 
Lil — Xi Ti4l41 — Liy1 
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En esta definición tenemos que asumir que 0/0 = 0. La curva B-spline está por 
tanto definida por una función spline polinomial que es de grado k y tiene derivadas 
continuas hasta orden k — 1. Además N;¿¡(t) > 0 y 


Observamos que la parametrización, al igual que en las curvas de Bézier, es poli- 
nomial y si k < n está definida a trozos. Hay curvas, como las cónicas, que no admiten 
una parametrización polinomial y por tanto no se pueden parametrizar mediante cur- 
vas B-spline, hablaremos de ellas en la Sección 3. Las curvas de Bézier son un caso 
particular de curvas B-spline, con grado n, el máximo posible, y vector de nudos igual 
a 

A A DE 


La curva B-spline está formada por trozos de curvas de Bézier, todas ellas del 
mismo grado k, el grado de la curva B-spline. El vector de nudos contiene el valor del 
parámetro t en los puntos de unión de los trozos de curvas de Bézier. 


Ejemplo 2.0.2. Calculamos con Maple, la parametrización de la curva B-spline de 
grado 2 con puntos de control By = (1,1), By, = (3,5),B2 = (5,3),B3 = (7,4) y vector 
de nudos X = (11,13, 13, T4,T5, 16,17) = (0,0,0,1,2,2,2). Es una parametrización 
polinomial a trozos 


P(t) =N ¡3(t)Bo + Nag(t)B1 + N3g(t)B2 + Nag(t) Ba 


SS (Eé+4t+1,-5e4+8t+1) si0<t<l, 
1. (3+4,8—6t + 2t?) E 


siendo N;z(t) las funciones de la gráfica de la figura 2.1. Observamos que la curva es 
la unión de dos curvas de Bézier, que se unen en el punto C' = (4,4). La curva de 
Bézier determinada por los puntos de control By, B,,C' tiene parametrización 


Pi(€) = (1— 0)Bp + 2t(1 —t)B, +90 = (- + 4t+1, 5 +8t+1), t € [0,1]. 


La curva de Bézier determinada por los puntos de control C, Ba, B3 tiene parametri- 
zación 


Pa(t) = (1— 1)Bpo + 2t(1 —t)B, +40 = (t? + 2t + 4, 21% — 2£+4), te [0, 1]. 
Si reparametrizamos para que t € [1,2] tenemos 


Pot - 1) = (3+1%,8-— 6t + 21%), + € [1,2]. 
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Figura 2.1: (Izqd) Funciones N,z(t), ¿1 = 1,...,4 de grado 2 de la base B-spline y 
(dcha) curva B-spline de grado 2, para cuatro vértices 


2.1. Vector de nudos 


El vector de nudos X', es un lista creciente de números reales, con multiplicidad a 
lo sumo el grado k + 1. La longitud de X depende del número de puntos de control 
n+1 y del orden k +1, siendo 


X 1 [21, ... Edel: 


Hemos comentado que los valores que aparecen en el vector X son valores de los 
parámetros en los puntos de unión de los trozos de curvas de Bézier que forman la 
curva B-spline. Cuando los nudos están uniformemente espaciados se dice que la curva 
es B-spline uniforme y en otro caso no uniforme. 


[1, 2,3, 4,5, 6] uniforme para n =2,k = 2, 
[1,3, 4,7,8, 10] no uniforme, 


[1,2, 2,3, 4, 4] periódico no uniforme. 


Los valores duplicados del nudo en la mitad de la lista hacen que una curva B-spline 
sea menos suave. Se dice que el vector de nudos es no uniforme periódico. En un caso 
extremo, un nudo de multiplicidad k+ 1 en la mitad de la lista de nudos significa que 
hay un lugar en la curva B-spline que se puede doblar en un punto de torsión. Por 
esta razón, a algunos diseñadores les gusta agregar y quitar nudos y luego ajustar los 
puntos de control para hacer curvas más suaves o figuras torsionadas. 

Si una lista de nudos se inicia con un nudo de multiplicidad k + 1, le siguen 
nudos simples y termina con un nudo de multiplicidad k +1 y los valores se espacian 
uniformemente, se dice que el vector de nudos es uniforme abierto (en inglés open 
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Figura 2.2: (Izq) Curvas B-spline de grados k = 1,2, 4,7. (Dcha) Código Grasshopper, 
haciendo visible el vector de nudos 


uniform), pero en general no tiene por qué serlo. 


[0, 0, 1, 2,3, 4, 4] uniforme abierto para n =4,k= 1, 

[0, 0,0, 1,2,3,3,3] uniforme abierto para n = 4,k = 2, 

[0, 0,0, 1/3,2/3, 1,1, 1] normalizado al intervalo [0, 1], 
[0,0,0,0,1,2,2,2,2] uniforme abierto para n=4,k=3, 
[0,0,0,1/2,1,1,1,1] normalizado al intervalo [0, 1], 
[0,0,0,0,0,1,1,1,1, 1] uniforme abierto para n =4,k=4. 


Se muestran en la figura 2.2 los vectores de nudos que toma Grasshopper para definir 
una curva B-spline de distintos grados. Observamos que en Grasshopper el vector de 
nudos tiene longitud n + k y los nudos aparecen con multiplicidad a lo sumo k. 


2.2. Influencia del grado 


En general, una curva B-spline de grado k con n+1 puntos de control está formada 
por n +1 k trozos de curva de Bézier de grado k. Fijado un polígono de control 
podemos obtener distintas curvas B-spline variando el grado. El grado mínimo es 
k = 1 y la curva B-spline es el polígono de control. Conforme aumenta el grado, la 
curva aproxima peor al polígono de control, (aunque aumenta la continuidad, véase la 
Sección 2.3). Observamos que modificar la longitud del vector de nudos es equivalente 
a cambiar el grado. llustramos estos fenómenos en la figura 2.2. 

Las curvas B-spline nos permiten tener control local. Si el grado k es menor o igual 
que n— 1, al mover un punto de control, el cambio no afecta a toda la curva B-spline 
sino a un número de vértices adyacentes. En la figura 2.3 se muestran dos curvas que 
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Figura 2.3: (Izq) Dos curvas B-spline de grado k = 2. (Dcha) Proyección al plano xy 
de ambas curvas 


tiene todos sus puntos de control iguales excepto el Ba, que en la curva plana es igual 
a (2,—1,0) y en la alabeada es (2, —1,2). Las curvas son de grado k = 2 y tienen 8 
puntos de control, por lo tanto están formadas por 6 curvas de Bézier B;,...,Bg de 
grado 2. Los 5 puntos de intersección Cy,...,Cs de la curva plana con su polígono 
de control determinan los puntos de unión de las curvas de Bézier. La curva B, tiene 
puntos de control C¡, B>,Ca3 y el desplazamiento del punto B» afecta también a las 
curvas adyacentes (afecta a Bi, Ba, B3). 

Las curvas B-spline admiten modo abierto y modo cerrado. Hasta ahora hemos 
utilizado el modo abierto, en el que la curva B-spline contiene los puntos inicial y 
final (Bo y B,,) del polígono de control. En el modo cerrado, el polígono de control es 
cerrado, Bo = B,, y la curva B-spline es una curva lisa que (en general) no contiene 
a dichos puntos. En el ejemplo de la figura 2.4, ilustramos que el componente Nurbs 
de Grasshopper permite elegir modo cerrado, con entrada C' igual a true. En modo 
cerrado tendríamos 9 puntos de control, ya que el primero se repite al final. En este 
caso hemos tomado grado 6. Observamos como la curva abierta tiene un pico en el 
punto By = By y sin embargo la curva cerrada es suave. 

Por estar formadas por una unión de curvas de Bézier, las curvas B-spline heredan 
propiedades de estas: 


1. Al igual que las curvas de Bézier, las curvas B-spline están contenidas en la 
envolvente convexa de los puntos de control. 


2. El ejemplo de la figura 2.3 ilustra que las curvas B-spline también son invariantes 
mediante transformaciones afines, en este caso observamos como la proyección 
de la curva de Bézier espacial sobre el plano xy, coincide con la curva B-spline 
de los puntos de control proyectados. 


3. La parametrización es de clase C*=1 (en los nudos). 
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Figura 2.4: (Izq) Dos curvas B-spline de grado k = 2, una abierta (en verde) y otra 
cerrada. (Dcha) Código Grasshopper 


2.3. Continuidad geométrica 


Sea T' la unión de dos (trozos de) curvas ', y P'3 en un punto extremo P común 
a ambas. Analicemos la continuidad de I' en el punto P. 


1. Por unirse en un punto extremo común P, el segmento I' es continuo en P, dire- 
mos que es C% (la continuidad clásica). También se dice que TP tiene continuidad 
geométrica G% en P. Las nociones de continuidad C% y G” coinciden. 


2. Si los vectores tangentes v¡ a I', y va a Pz en P tienen la misma dirección v, = 
aw, a > 0, diremos que T' tiene continuidad geométrica G* en P. Observamos 
que si la parametrización P(t) de T' es de clase C* en P entonces los vectores 
tangentes tienen la misma dirección y la misma magnitud, por lo que C! es más 
exigente que G*. Así, C* implica G* pero no al revés. 


3. Teniendo continuidad geométrica G* en P, diremos que T' tiene continuidad 
geométrica G? si los vectores curvatura k1 y ka en P de I'; y Ta tienen la misma 
dirección k¡ = QGk2, Y > 0. 


Aunque en muchas aplicaciones es suficiente exigir continuidad geométrica G”, existen 
otras (aquellas que dependen de la suavidad de las curvas o superficies) en las que 
es necesaria una continuidad G? o incluso C?, especialmente las que dependen de 
una transición suave de la luz reflejada, como por son por ejemplo los automóviles 
[Rogers]. 


Ejemplo 2.3.1. Conectamos dos curvas l', y > de Bézier mediante el componente 
Connect de Grasshopper, que toma como input una lista de curvas C =4T,,TP2) y un 
número entero G, la continuidad geométrica 0,1 6 2. Las curvas has sido creadas con 
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Figura 2.5: Curvas conectadas con continuidad (arriba izq) G%, (arriba dcha) G?, 
(abajo dcha) G?. (Abajo izq) código Grasshopper 


el componente Nurbs, la primera de grado 2 (con 3 puntos de control) y la segunda de 
grado 3 (con 4 puntos de control). Utilizamos el componente CruGraph para visualizar 
el gráfico de curvatura de la curva resultado T. El componente Ins nos permite in- 
sertar las curvas Py y Pa en una misma lista. En la figura 2.5 mostramos el resultado 
de conectar las curvas utilizando los distintos tipos de continuidad geométrica G%, G* 
y G?. 

Es interesante observar que, para conectar T', y T'2 con continuidad G*, Grasshop- 
per genera una curva B-spline de grado 3 con vector de nudos periódico 


x=/(00,0,1,11.32233,3] 


Si le pedimos continuidad G?, genera una curva B-spline de grado 5 con vector de 
nudos 
X =[0,0,0,0,0,1,1,1,1,1,2,2,2,2,2,3,3,3,3,3]. 
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3 Curvas B-spline racionales (NURBS) 


Ya hemos hablado de casi todos los ingredientes de las curvas NURBS (Non uni- 
form rational B-spline), nos falta explicar qué significa ser racional. Es preciso matizar 
que el nombre NURBS no es excluyente, incluye a todos los casos anteriormente estu- 
diados. Son curvas NURBS las curvas de Bézier y las curvas B-spline, tanto uniformes 
como no uniformes. El aspecto racional es el que nos permitirá dibujar todas las cóni- 
cas mediante este tipo de curvas, véanse los detalles de la Sección 3.1. 

Una curva NURBS plana es la proyección central de una curva B-spline espacial 
sobre un plano. Análogamente una curva NURBS espacial es la proyección central, 
sobre un hiperplano, de una curva B-spline que vive en un espacio de dimensión 4. 
De forma más precisa, sean Bo,..., B,, puntos de control en R%, supongamos que 
B;= (a;,b;,c;). Tomemos coordenadas homogéneas B" = w;(1, as, b;,c;) ER*, w, > 0 
de B;. Una curva B-spline racional espacial es la curva afín de una curva B-spline 
proyectiva 


n+1 


PO) = Y BR Nit), a <t<b. 
i=l 


Pasando a coordenadas afines obtenemos la parametrización de la curva afín, la curva 
B-spline racional 


P(t) EE A wi 1B; IN +1 (0) 
E wi-1N ¿41(8) 


ALEC (3.1) 


Llamamos peso del punto de control B; a la primera coordenada w, de BP. Nos 
restringimos a pesos positivos w, > O para evitar (en este manuscrito) utilizar puntos 
en el infinito. El incremento del valor del peso w; acerca la curva NURBS al punto 
de control B;, si el valor del peso decrece se produce el efecto contrario y la curva 
NUREBS se aleja del punto de control. Este desplazamiento tiene un efecto local en la 
curva NURBS. 

Las curvas NURBS que tienen todos sus pesos iguales, w, =w, ¿+= 0,...,n son 
curvas B-spline, la parametrización (3.1) es polinomial ya que el denominador es 
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constante 


n+1 n+1 
y wi-1N¿k(t) =u y Nir(t) =w, 
i=1 ¿=1 
debido a que la suma de las funciones B-spline NV; y(t), += 1,...,n+1 es igual a 1. En 


este caso, la curva parametrizada por P(t) sería una curva B-spline en el hiperplano 
w =w, semejante a su proyección en el hiperplano w = 1, parametrizada por P(t). 
En otro caso, si existe algún peso distinto de los demás, la parametrización P(t) es 
racional y tenemos una NURBS genuina, que no es B-spline. 

El diseñador tiene que tener en cuenta que para modelar utilizando curvas NURBS, 
necesita establecer el valor de: 


1. Los puntos de control Bo,..., B,. 
2. El grado k menor o igual que n y el vector de nudos X. 
3. Los pesos w; > 0,2=0,...,n. 


Podría jugar con todos estos datos hasta obtener la curva deseada. Si utiliza curvas 
B-spline, sólo podría modificar puntos de control, grado y vector de nudos, los pesos 
son todos iguales a 1. El caso menos versátil es el de las curvas de Bézier, para las 
que sólo podría modificar los puntos de control, el grado es fijo (número de puntos de 
control menos uno, y el vector de nudos también) y los pesos son todos iguales a 1. 

El componente de Grasshopper que nos permite controlar todos los ingredientes 
de una curva NURBS es NurbCrv, que toma como entradas: 


1. P la lista de puntos de control Bo,..., B,. 
2. W la lista de pesos wo,...,Wp. 
3. K el vector de nudos X. 


El grado de la curva NURBS queda determinado por la longitud 1 del vector de 
nudos, que en Grasshopper es l =n +k (aunque en la definición de la base B-spline 
es n + k +2). Como ya comentamos en la Sección 2.1, modificar el vector de nudos 
equivale a elegir el grado. 

En el ejemplo de la figura 3.1 ilustramos el efecto de cambiar el valor del peso 
w1. La curva de partida es una circunferencia del plano xy, de radio 20 y centrada 
en el origen. Observamos que es una curva B-spline racional, formada por 4 arcos de 
curva de grado 2 (cada arco sería una curva de Bézier racional). El arco del primer 
cuadrante está determinado por los puntos Bo, B, y Ba. 
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NURBS dadas por: puntos, pesos y nudos 


'Camblar pesos para acercarnos a un nodo Ñ 


Figura 3.1: (Izq) Una circunferencia y la curva NURBS que resulta al cambiar el peso 
de un punto de control. (Dcha) Código Grasshopper para modificar puntos, pesos y 
nudos en una curva NURBS, en el componente NurbCrv 


3.1. Cónicas 


Las secciones cónicas son curvas descritas por ecuaciones de segundo grado, que 
además tienen parametrizaciones racionales. Por este motivo es interesante analizar 
las curvas NURBS dadas por 3 puntos de control By, B,, Ba y con grado máximo k = 
2. Veremos que son arcos de cónica, y que sólo las parábolas tienen parametrizaciones 
polinomiales. 

Supongamos que trabajamos con puntos en el plano B, = (a;,b;,), ¿ = 0,1,2. 
Elegimos coordenadas homogéneas de dichos puntos 


B' = (w;, w,0;, 0:b,), 1 = 0, 1,2, 


con pesos w, > 0. Los puntos Bf, Bf y Bl determinan una curva de Bézier C de 
grado 2 espacial, parametrizada por 


P() =» Jas(t)Bi = 


O A A 
(Polt), PL(t), Palt)). 


La curva C es una cónica proyectiva, cuya cónica afín C está parametrizada por 


(PLE Patty _ a wiJa:(t)B, 
AU o 200) M pd wit: (t) 
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Figura 3.2: Curva proyectiva y su proyección sobre el plano afín (xy = 1) 


Ejemplo 3.1.1. Tomemos puntos de control By = (0,1), B¡ = (1,1), B2= (1,0) y 
pesos Wy = 3, w = 5, wi = z. La parametrización de la cónica proyectiva (curva de 


Bézier de grado 2) es 


5 1.11 2 1 3 
Pi =1| -=* t == E 
(+) (5 MEJESE: 3 +3 ra 


Representamos en la figura 3.2 la cónica proyectiva C y su proyección central C sobre 
el plano afín (si denotamos por (xp, 11,12) las coordenadas de un punto del plano 
proyectivo, estamos proyectando sobre el plano ry =1). 


Observamos que la primera coordenada de P(t) es Pol[t) = N_, Ja :(t)w, y los 
puntos del infinito de C son aquellos para los que Po(t) = 0. Es conveniente tomar 


W=wW=l y w =u. 


Estudiemos cuántos puntos del infinito tiene C en función del valor de w, lo que 
equivale a determinar el tipo de cónica. Las soluciones de 


Polt) = (1—t)?+2(1 — t)tw+t?=0. 


son 


_ —l1+w+y-1+u? 


t, = 


o Li my 13 ue 
2(w 1) da Mw=1D) 


Fácilmente observamos que la cónica proyectiva C tiene los siguientes puntos en el 
infinito: 


pl 


1. Siw= 1 entonces t¡ = ta, un punto doble en el infinito. La cónica afín C es una 
parábola. 


2. S10 < w< 1 entonces t; y tz son valores complejos, dos puntos del infinito 
complejos. La cónica afín C es una elipse. 


3. Si w > 1 entonces t, y ta son valores reales, dos puntos del infinito reales. La 
cónica afín C es una hipérbola. 


Podemos calcular la ecuación implícita f(x, y) = 0 de la cónica afín eliminando la 
variable t de los polinomios 


Polt)jz — Pi(t) =0, 
Polt)y = Pat) = 0. 


Para ello podemos utilizar el comando Maple, resultant. Por ejemplo, si tomamos 
Bo = (0, 1, B; = (1,1) y Bo) = (1,0) 


Pi(t)= 2wt =2w0 +1 =2+ 2, 
Polt) = 1? + 2wt — 2wt?. 
y así 


Fx, y) =resultant(Po(t)x — Py (t), Polt)y — Pa[t), t) 


=1 — 2y — 21 — 4w* + y? + 2? 4 2yx + 4xw? — 4yw*z + 4yw? = 
1—4u? —14+24? —1+2u? 1 

=(1 x y)| -1+2? 1 1 — 2? a 
ALI LD 1 y 


Observamos que la matriz de la cónica es degenerada sólo cuando w = 0. Además si 
w= Ya entonces C es el arco de la circunferencia 2? + y? = 1 con parametrización 


p(t) = 2? +tyV2-y2 tv2-42+1-2+t 
22 + 42-122 +1-2t' 22 + t/2-1t2y/2+1-— 2t 


Ejemplo 3.1.2. Tomemos puntos de control By = (0,0), B, = (3,4), Ba = (4,0) y 


pesos Wy = 1, w, = w, wa» = 1. Representamos en la figura 3.3 tres arcos de cónica 
1 


para distintos valores de w. Siw = es un arco de elipse, w =1 un arco de parábola 


y w=3 un arco de hipérbola. Las parametrizaciones son: 
t1+3t  4t—4t? 1 
(21) A 
(-28* + 6t, 8t — 81?) para w=l, 
141? + 18  24t — 24t? 
(+ 1 42 + 4t4 1 


) para w= 3. 
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Figura 3.3: Cónicas para (de abajo a arriba) w = » 


Nótese que sólo la parábola tiene parametrización polinomial. 
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